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Abstract 

Fermions localized within vortex cores can form one-dimensional 
Fermi liquids. The nonzero density of states in these Fermi-liquids can 
lead to instability of the symmetric structure of the vortex core. We 
consider a symmetry breaking which is obtained due to spontaneous 
admixture of the spin-triplet p-wave component of the order parameter 
in conventional s-wave vortex in the presence of magnetic field. This 
occurs at low enough temperature and leads to an asymmetric shape 
of the vortex core. Similar phenomenon of the symmetry breaking 
induced by the core fermions occurs in electroweak Z-strings ||. 

Introduction. Recently new techniques have been developed, that in 
principle allow to probe the core structure of the individual Abrikosov vor- 
tex in superconductors. This includes a scanning tunneling microscopy [[[]], 
electron holography, Lorentz microscopy, etc. (see @ and Refs. there). That 
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is why the theoretical investigation of the symmetry of the core and of the 
energy spectrum of the fermions localized in the core is now up-to-date. 

Three possible types of the energy spectrum E(Q, k z ) of the fermions 
localized within the vortex core have been discussed in Ref . || . (i) The spec- 
trum E(Q, k z ), where Q is the generalized angular momentum, has a gap on 
order A 2 /Ep Q (here A is the superconducting gap and Ep is the Fermi en- 
ergy in the normal Fermi liquid), (ii) One or more branches of the spectrum 
with particular Q's cross the zero level as functions of the linear momentum 
k z at some points k z = kp^Q. In this case fermions occupying negative energy 
levels form one-dimensional (ID) Fermi liquids with the Fermi points at kF,Q 
(see also an example of the continuous vortex in 3 He-A [|5]. (iii) The flat 
band where all fermions have exactly zero energy in the finite region of the 
momenta k z , which is called the Fermi condensate ||. 

It was emphasized that both the ID Fermi liquid and ID Fermi condensate 
can be unstable towards a further breaking of the vortex symmetry at low 
temperatures, in particular the vortex can become nonaxisymmetric or the 
vortex line can transform into a spiral |J (see also [[?[]). 

Recently such kind of instability of the vortex caused by the fermion 
zero modes in the core has been discussed for the Z-string solution in the 
Weinberg-Salam model of electroweak interactions ||. In this case the sit- 
uation corresponds to the case of ID Fermi liquid. While in the Ref.0 it 
was stated that this can lead to the absolute instability of the Z-string, the 
situation is actually not dramatic for the existence of the locally stable Z- 
string and is similar to that which has first been discussed by Peierls [[J. In 
the Peierls model the ID electron liquid coupled to the lattice is unstable 
at low temperature towards the dimerization of the lattice. In the Z-string 
the ID massless relativistic fermions coupled to the bosonic field of the order 
parameter (Higgs field) lead below some critical temperature to the sponta- 
neous symmetry breaking, which is realized, e.g., in the development of the 
small amount of the upper component of the Higgs field or of another mode 
of instability. The transition temperature T c core of such instability is usually 
exponentially small. 

We discuss here the similar instability for the condensed matter vortices. 
We have found that in most cases the transition temperature T c core is ex- 
tremely small of the order Aexp(— Ep/A), but in some cases T c core can be 
reasonable, of the order A 2 /Ep. 

s-wave superfiuids. Let us consider as an example the situation in 
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superconductors with s-pairing. We suppose that the particles are chargeless, 
and therefore we neglect all magnetic effects. In the case the dependence of 
the order parameter on coordinates around the axis of one-quantum vortex 
is given by A(r) • exp(i0) in cylindrical coordinates. The magnitude of the 
gap A(r) equals zero at the axis r = and tends to the bulk value A for 
large distances r ^> £, where £ = Vp/A is the coherence length. 

The system is invariant under translations along the axis z. If the spin- 
orbit interaction is neglected there is also the SU(2)s group of spin rotation, 
which is reduced to 5*0(2)5 subgroup if the magnetic field is applied. The 
"axial" symmetry of the order parameter in the vortex is characterized by 
another SO(2)q with a generator which is the combination of the orbital 
momentum L z and the particle number operator I [[TTJ: 



where m — 1 is the winding number of the vortex. The fermion eigenstates 
are characterized by the quantum numbers k z , Q and the direction of the 
spin (| or |). The energy spectrum of low-lying excitations is given by [§]: 



with Uq ~ A 2 / Ep for k z <C kp and u)q — > oo for k z —>■ kp (this approximation 
is true for uj ^ A) @|. For m = 1-quantum vortex Q can assume only 
half-integer values (L z is integer, and I is +1 and —1 for particles and holes 
correspondingly) . 

Let us apply a magnetic field H parallel to the axis of the vortex. Then 
the energy spectrum is shifted down (up) for fermions with spin up (down): 



and we get the picture (fig.|T|) where the branches EQ^{k z ) and EQ t i(k z ) 
intersect for jj,H > Qojq{Q)/2. 

The fermions near the points kp^Q can be considered as a one-dimensional 
Fermi-liquid, and kp^q plays the part of the Fermi- momentum. The nonzero 
density of states can lead to a new pairing of the one- dimensional fermions, 
the character of which depends on the sign and the magnitude of the interac- 
tion of ID fermions in different channels. Let us show that both symmetries 




(1) 



E(Q,k z )=Q-u {k z ) 



(2) 



E hl (Q, k z )=Q ■ u (k z ) =f jiH, 



(3) 
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SO(2)s and SO(2)q will be broken at low enough temperatures, if the in- 
teraction in the spin-triplet p-wave channel is attractive, and the core of the 
vortex can become anisotropic. 

p-wave component in the core. In our model an extra interaction of the 
fermions in ID Fermi liquid results from the interaction of the initial bare 
fermions in the p-channel and we are looking for the instability of the vortex 
core towards nucleation of the admixture of the spin triplet p-wave component 
in the core of the type: 

5A = g A(d • <?)i<3- 2 (a • k) f(r) exp(iiV0) . (4) 

As in Ref.[|8| we assume that f(r) is concentrated in the core since it should 
influence only the core fermions. f(r) is of the order of unity in the core 
and decreases to zero outside the core of the vortex (i.e., at r ~ £); in 
addition /(0) = 0, if iV 7^ 0. The parameter A is the bulk value of the 
s-wave gap and the dimensionless parameter g is thus the magnitude of the 
admixture of the p-wave order parameter in the core region as compared to 
the s-wave component. For simplicity we can assume a trial function of the 
form f(r) = exp(— r/£) for N = 0. The spin vector d should be orthogonal 
to the magnetic field to break the spin rotation symmetry. 

The suggested change of the order parameter breaks Q-symmetry as well. 
To show this note that 5 A is the sum of three orbital harmonics: a ■ k = 
a + /c_ + a-k + + aok z where k± = k x ± ik y . We suppose that only one a 
is nonzero, namely a\ with some value I = 0, ±1 of the ^-projection of the 
Cooper pair orbital momentum. It is easy to check that new order parameter 
is not Q-invariant for N + I ^ 1 and this means the violation of the axial 
symmetry of the whole core structure, g is the dimensionless parameter 
which shows the strength of the deformation. This symmetry breaking leads 
to mixing of Q, j and —Q, j fermions states in the same way as u and d quarks 
are mixed in the Z-string in the Naculich scenario of the instability of Z- 
string. Due to the noncrossing theorem for the levels of the same symmetry, 
the two branches which crossed each other at kp,Q repel each other and the 
gap between the states appears which is proportional to g. 

Deformation of the fermionic spectrum. Without the admixture the 
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Bogolyubov-Nambu hamiltonian is given by 4 x 4 matrix 
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( t 
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-£ y 



(5) 



Here £ = (l/2m)(—d 2 /dr 2 — kp). The eigenstates of Ho with energies (0) for 
spin up and distances r 3> fc^ 1 are given by M 



Xk z ,Q,-[ 



( XI \ 




v X2 y 



amst-e ik ' z e iQ *H v (qr)e- K ( r) ' 



^ exp(|(0 + V(r))) ^ 



V -?exp (-§(</> + ,/ 

(6) 



where ifj, is the Hankel function with index z/ = yQ 2 + 1/4; q 2 = kp — k 

r 

Kir) = — [ A(r')dr', 
a J 



exp(2(K(r) - K(r'))) ( — + Q 



V 9 
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For spin down the eigenstates are 



Xk z ,Q,l 



( \ 

Xi 

~X2 

V o y 



(7) 



The change 5 A leads to the perturbation of the mean field Bogolyubov- 
Nambu hamiltonian. The angular integral in the matrix element of this 
perturbation between the states Xk z ,Q,i and Xk z ,-Q,i does not vanish only for 



N + l + l±2Q = 



(8) 



(i.e. iV + I should be an even integer). In this case the matrix element 
between the spin- up and spin-down states is of the order of gA. For other Q 
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the integral over d<\> in the matrix element vanishes. The 2x2 hamiltonian 
for two mixing states (k z , Q, |) and (k z , —Q, j) for k z near k F Q is then given 
by: 

( E old gA \ 
{ gA —E id I ' 



(9) 



Here E Q i d = v F; Q(k z — k F ^), the derivative vp,Q = Qoj' q of the energy over 
momentum k z represents the "Fermi" velocity of the ID Fermi-liquid. The 
modified energy spectrum takes the form 



E new = ±^v% Q {k z - k F , Q y + g 2 A 2 . (10) 

At zero temperature fermions occupying the states with negative energy near 
kF,Q gain in energy after the change in the order parameter. We are going 
now to compare energy gain and energy losses due to the perturbation g. 

Energy balance and zero-temperature deformation of the vortex core. The 
main contribution to the difference between the energy of the fermionic vac- 
uum in symmetric and nonsymmetric states comes from the logarithmic term 
which arises due to appearance of the energy gap in ID Fermi-liquid: 

•I-* gain J W~ \ new -^old) 

S 2 A 2 k fd{k-k F , Q ) 



v f,q J \k-k F)Q \ 
fci 

The lower cut-off is determined by the gap itself 

fci = — . (12) 

V F,Q 

The most interesting situation takes place for k F> Q not too close to k F , in 
this case the Fermi- velocity v f> q is small and the energy gain increases. This 
occurs in field H near its threshold values H C (Q) at which the branches Q, j 
and —Q, j start to cross each other and the "Fermi- momentum" k Fy Q appears 
for the first time. The "Fermi-momentum" k Fy Q is given by the condition 
Quj (k F) q) = fiH, that is we are interested in the fields fiH ~ QA 2 /E F . In 
this region the upper cut-off in the logarithmically divergent integral in (O) 
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is &2 ~ kp,Q ~ kpJ Qu H (o) ~~ This l ea ds to the result 



M)(o) 

On the other hand the energy loss due to the perturbation of the su- 
perconducting order parameter in the core of the vortex is of the order of 
Ei oss = Epkp g 2 /ji- Here 71 is the relative magnitude of the attractive in- 
teraction of bare fermions in the p-channel compared to that in s-channel. 
For N ^ the contribution to the gradient energy of space inhomogeneity of 
SA increases by the amount ~ EpkpN 2 g 2 . This consideration is valid only 
if the latter term is smaller than the former one. 

Comparison of the energies shows the instability, the relative magnitude 
g of the new p-wave order parameter at zero temperature being of the order 

9 ,A IQI (y ! exp (_ IQI ^ ( I + Jv2) ). (14) 

Here we omitted the coefficients of the order of unity. N(Q) is given by (|8|). 
At small Q the exponent exp(— 0(1)) is not too small for kp,Q ~ kp, and 
the transition temperature is reasonable. On the other hand the numerical 
factor in the exponent 0(1) can be large. In this case (even for large Q) we 
can take kp^/kp to be small which leads to the exponential increase in the 
magnitude of g followed by just power-law decrease in the preexponential 
factor. 

Symmetry breaking scheme. The above result means that there should 
be a second-order phase transition at T c core into the broken symmetry state. 
This is the state with the additional superfluid order parameter (01) with 
p-pairing. The k-dependence of this order parameter is described by the 
quantum number / = 0, ±1, and the spatial ^-dependence by another integer 
N = ±20, — / — 1. So, there are six competing structures (for three values 
of I) with comparable transition temperatures. They correspond to different 
symmetry groups of the new vortex core structure below T c core . 

The spin structure of the new order parameter is described by the vector 
d in the xy-plane. For d = x±iy the symmetry SO(2)q x SO(2) s is reduced 
to the combined symmetry group 5 , 0(2)q± 2 |q±i|s- For other d in the plane 
the symmetry is reduced to the discrete group Z±\q±i\ the elements of which 
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are Q-rotations accompanied by the possible inversion of d. In the former 
case, say, density is still axially symmetric, and in the latter case it is not. 

Transition temperatures. The temperature T c core at which the transition 
to the asymmetric core state occurs is 

T c core ~ 9 A . (15) 

(see the phase diagram, fig.@). The maximal transition temperature in the 
Q-th region and the range of the fields where the bare structure is unstable 
(the height and the width of the region in H — T-plane) decrease with Q: 

T c , max (Q) ~ A»H(Q) ~ ^ig. (16) 

The minimal field is of order (A 2 /Ep)/ fi. This field can be less than the 
upper critical field H C 2 if the effective fermion mass m* is larger than the 
bare electron mass m. 

Vortex in 3 He-B. In the spin-triplet superfluid 3 He-B the situation is al- 
most equivalent to that in s-superfluids. In the most symmetric B-phase vor- 
tex the branches of the energy spectrum are described by quantum numbers 
k z , Q and helicity A = ±1 @. The branches _Eq =0A=1 (/c 2 ) and EQ =0X =-i(k z ) 
intersect at the the point k z = 0, E = even in zero external magnetic field. 
The branches Eq X =i and E_q^=-i begin to intersect at k z ~ kp in the fields 
fiH c (Q) ~ Qcu ~ QA 2 /E F (fig.|). 

In the case there is no conservation of longitudinal spin component S z , 
and the perturbing order parameter may have vector d of spin anisotropy 
parallel to z-axis. The deformation of the vortex can occur within the triplet 
pairing state; the singlet pairing (e.g., s-pairing) is accepted as well. 

Transition temperatures, critical fields and the whole phase diagram are 
more or less the same as for s-superfluids. For the two branches mentioned 
above which intersect in zero field the transition temperature is of the order 
of T c ~ A 2 / ' Ep. It was found experimentally and theoretically that the most 
symmetric vortex in 3 He-B is unstable towards the symmetry breaking at 
least at high temperature of order T c JTUL [□], O. Thus if the symmetry is 
restored at lower temperature it is broken again at very low T. 

The states with definite A are the mixtures of states with different particle 
number and different S z = ±1/2. The analysis shows that in the Q-th region 
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the angle dependence of the perturbation of the order parameter can be given 
by an integer 

N = ±2Q -l-m + E 

where S = 0, ±2. Therefore for each I we get three different values of N and 
thus several competing types of the symmetry breaking. 

Conclusion. The vortex structure, which contains the one-dimensional 
Fermi-liquid of the core fermions, is unstable towards the symmetry breaking. 
In the case of Abrikosov vortex in the s-wave superconductors this leads to 
the breaking of the axial symmetry of the vortex core in some regions in 
H — T-plane. 

Similar symmetry breaking occurs in the electroweak Z-string, where the 
mixing of the d and u quarks opens the gap in the fermionic spectrum ||. In 
the Abrikosov and the 3 He-B vortices this corresponds to the hybridization 
of Q, | and —Q, j states of fermions. For the Abrikosov vortex this leads to 
appearance of the p-wave pairing component in the core. 
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Figure 1: Spectrum of fermions in the vortex in s-wave superconductor. The 
levels are shifted by magnetic field. Some branches cross the zero energy level 
forming ID Fermi liquid. The fermions with opposite spins on the Fermi level 
are paired at low T which leads to broken symmetry of the core. 



Figure 2: Schematic dependence of the temperature of the core transtion on 
the magnetic field. Within the shaded area the symmetry of the vortex core 
is broken. 



Figure 3: Spectrum of fermions in the most symmetric vortex in 3 He-B [[|. 
A = ±1 is helicity. 
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